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PREFACE 


This report presents the description of a computer code which used the 
spectrum characterization methods developed in a research program entitled 
"Improved Methods for Predicting Spectrum Loading Effects.'' This program 
was administrated by the Flight Dynamics Laboratory of the Air Force Wright 
Aeronautical Laboratory, Wright-Patterson Air Force Base, Ohio, under Contract 
F33615-077-C-3121, Project 2401, "Structural Mechanics" Task 240101, "Struc- 
tural Integrity for Military Aerospace Vehicles," Work Unit 24010120. R.M. 
Engle (AFWAL/FIBEC) was the Air Force project engineer. 


This research program was primarily conducted by personnel from the 
Fatigue and Fracture Mechanics Group, Dynamics Technology, Structure Systems, 
supervised by George E. Fitch, Jr., supervisor, Joseph S. Rosenthal, manager, 
and Dr. Leslie M. Lackman, director. James B. Chang was the program manager 
and principal investigator. Edward Klein participated in the original 
development of this computer code. The spectrum characterization procedure 
used in this computer code was developed in Phase I of this research effort 
by Drs. Masanobu Shinozuka and Rimas Vaicaitis. 
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Section I 


? INTRODUCTION 


Current military standard MIL-STD-1530A "Aircraft Structural Integrity 
Program, Airplane Requirements" 1) states that two major activities designed 
to focus attention on each potential crack problem shall be included in the 
force management task. They are the Force Structural Maintenance (FSM) 

Plan and the Individual Airplane Tracking (IAT) Program. In addition, force 
management activities also include the loads/environmental spectral survey 
(L/ESS), updating the design analysis, developing inspection and repair 
criteria, and forming a structural strength survey. 


Mate? 


The objective of the IAT program is to predict the potential crack growth 
in critical areas of each airframe, keyed to crack-growth limits, damage- 
tolerance limits, inspection times, and economic repair times. In the JAT 
program, an individual airplane tracking analysis method which establishes and 
adjusts inspection and repair intervals for each critical area of the airframe, 
based on the individual airplane usage data, must be developed suited for a 
particular aircraft system. The damage-tolerance and durability analysis and 
associated test data are used to establish the tracking analysis method. This 
tracking analysis provides the capability to predict crack-growth rates, time 
to reach the critical crack sizes, and crack size as a function of total flight 

> time and aircraft usage data. : 


According to a survey conducted by the University of Davton/Lockheed/ 
Vought team, there are about 11 IAT methods being used in the IAT programs of 
25 aircraft svstems in the Air Force inventory. Among them, five methods are 
vased on crack-growth analysis (2). The common practice is to employ 2 cycle- 
by-cycle crack-growth computer program to compute the crack growth for a 
flight or a number of flights. 


From the economical point of view, the use of a cycle-by-cycle crack- 
growth computer code to compute the crack growth in the IAT program is 
definitely not cost effective. This is not only because the cycle-by-cycle 
crack-growth analysis consumes too much computation time, resulting in 
excessive computer cost, but also because there is no need for an accurate 
representation of the crack-growth behavior on the stress-level-by-stress- 
level basis for performing individual airplane tracking. Furthermore, 
it is highly desirable to operate the crack-growth analysis code on a 
variety of computer systems. It is even more desirable that such codes can 
be operated on the onboard type of minicomputers. Capacities of such 
computer systems prevent the use of the sophisticated cycle-by-cycle growth 
method. 


i 


It is for those reasons a fatigue crack-growth prediction procedure 
which uses two random spectrum characterization methods developed in phase I 
of this research effort, "Improved Methods for Predicting Spectrum Loading 
Effects" © , was formulated. These two methods were identified as Method 
I and Method II in Reference 3. Method I is used to develop equivalent load 
spectra in terms of constant-amplitude stress histories (one cycle per flight). 
Method IT is used to characterize each mission segment in a flight by a 
constant-amplitude load segment; a flight is assumed to consist of constant- 
amplitude mission segments. 


A computer code, FLTGRO, was developed in the second phase of the 
aforementioned research effort. Both Method I and Method II have been 
implemented into FLTGRO. It has been subsequently used in providing analytical 
analytical predictions on 41 random flight spectrum test cases. Analytical 
predictions were correlated with the test data. Very good correlations were 
shown. Results of the 5st data correlations were documented in the final 
report of this program’ °. This report presents the description of the two 
spectrum characterization methods implemented in FLTGRO and provides detailed 
instructions to the user for executing the FLTGRO program. 


to 


Section I] 


TECHNICAL DISCUSSTON 


The crack-growth prediction methodolopy implemented in this computer code 
uses two spectrum characterization approaches, One approach, identified as 
Method I, is able to characterize the random flight spectrum or the ordered 
flight spectrum into the equivalent constant-amplitude loadings in one- 
evcle-per-flight format. The second approach, identified as Method I], is 
able to represent each flight segment in a flight such as climb segment, 
maneuver segment, descent segment, etc, by a constant-amplitude load 
segment. A flight is then assumed to consist of several constant-amp] itude 
segments as shown schematically in Figure 1. The following are brief des- 
criptions otf the necessary steps required for these two approaches. Refer 
ro“teference 1 for a detailed description of each method. 


2.1. METHOD [T-ONE-CYCLE-PER-FLIGHT CRACK-GROWTH PREDICTION METHOD 


a. Generate sample stress histories for all missions and mission 
segments considered in certain combinations. Repeat the procedure 
NA times to produce Na consecutive sample flight spectra which 
is called a unitblock. 


b. Use a cycle-by-cyele crack-growth analysis computer routine to 
evaluate the crack growth sa due to the unitblock flight spectra 
under a prescribed number of different values of initial crack 
$size a. Then, calculate da/dF = sa/Na. 


c. Establish the following relationship on the basis of the foregoing 
nunerical results: 


dajdl = ciky' 


where K is a measure of the stress intensity factor representing 
the overall effect of the unitblock on the crack growth. Ina 
mathematical torm, k is written as 


aby db 


Koes (a) 


3 


STRESS MANEUVER 
SEGMENT 


CLIMB 
SEGMENT 


DESCENT 
SEGMENT 


(a) RANDOM FLIGHT 


—— — TRANSITION CYCLE 


STRESS o (KSI) 


. FLIGHT —————————_— 


(b) MULTISEGMENT LOAD 


Figure 1. Spectrum Schematic 
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where (4.7 ee represents the statistical average of the bth power of 
the stress ranve ao in the stress history and $a) is a function of 
crack size and other geometries; e@.g., for ua center-through crack in 
a plate under tensile load, 


y(a) = Wea Vsec (sa/W) 


Plot da/dF against K for a number of different values of crack size 
aon a double logarithm scale. Then determine the two parameters C 
and i. Each of the flight spectra representing a particular combina- 
tion of the stress parameters with the rate of crack-prowth evaluated 
is then replaced by a constant-amplitude stress spectnm with the 
equivalent rate of crack growth. 


Repeat the same procedures for other combinations of stress param- 
eters, resulting in corresponding values of C, +, and (.-b)1/b and 
also in corresponding equivalent constant-amplitude stress spectra. 
Thus, obtain a number of sets of parameters 
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and a number of growth rates applicable to each of the particular 
stress parameter combinations. 
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where \N is the total number of stress-parameter combinations toa be 
considered. Hence, the problem of mission mix and mission sequence 
can be solved by choosing those growth-rate equations pertinent to 
the missions in the mix and applying them in the sequence correspond 
ing to the mix. 


The crack-growth life (F) in flights is obtained by the modified 
moditied linear approximation damage accumudation technique. 


2.2. METHOD II - MULTISEGYENT-PER-FLIGHT CRACK-GROWIH METHOD 


The one-cycle-per-flight approach may not be appropriate for the highly 
maneuverable fighter-type aircraft. Hence, a multisegment-per-flight crack 
growth methodology was developed. Figure 1 schematically illustrates the 
conversion of the random flight spectrum to the multisegment-per-flight 
spectrum, 


The multisegment-per-flight crack-growth methodology also uses the 
analytical procedure as does the one-cycle-per-flight, except the crack- 
growth rate is on a cycle-per-segment (da/dS) basis rather than a cvcle-rer- 
flight (da/dF) basis. The following are brief descriptions of the necessary 
required steps: 


a. Generate sample stress histories for a specific flight segment such 
as a Manuver segment; repeat the procedure “a times to produce *y 
consecutive sample flight segments (a uniblock). 


Use a cycle-by-cycle crack-growth computer program to calculate the 
crack growth 4a due to the unitblock flight segment starting from a 
prescribed initial crack size a9; then, determine da/ds = Aa/n;. 
Repeat the procedure for a predescribed number of different ao 
values. 


Establish the following relationship on the basis of the aforemen- 
tioned results: 


da/ds = c¢ (ky 


where kK is the measure of the stress intensity factor representing 
the overall effect of the unitblock on the crack growth. Ina 
mathematical form, K can be represented by 


Kee” ¥(a) 
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The term Gone? is the statistical average of the b-th power 


of the stress range, Ac; ¥(a) is a function of crack size and 
other geometrica] parameters. 
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d. Plot da/ds against kK on a double logarithm scale; then, determine the 
two parameters C and +. Choose either the appropriate value of 
q omax OY cmin and then calculate the corresponding smin OY “max Pascd 
on (.c5)1/b value, resulting in the corresponding value of C, +, and 
Smax,s “min of the equivalent constant-amplitude flight segment. 


e. Repeat the same procedures for other tlight segments and thus 
obtain a number of sets of parameters 


and a number of growth rates applicable to each of the particular 
flight sepments 


\ 
BS alse. Aoee ae ck 
(da/ds} Cc) (kK, ) 
pan A, 
(da/ds), = C, (K,) ~ 
pee hs | 
(da/ds) ; 7 Ce (k;) 


where i is the number of flight segments considered in a typical a 
flight. 


f. Choose the number of cycles nj in each segment and then convert the 
crack-growth-per-segment rate (da/ds) into the crack-prowth-per-evele 
rate (da/dn), resulting in 


(da/dn) ; = G-| (da/ds) = Gr 
1 1 


g. Use a cycle-by-cycle crack-growth analysis computer program to 
calculate the growth behavior. 
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The evele-by-cvele crack-growth subroutine built in FLTGRO is identified 
as ESTCAL. The crack-growth methodology implemented in ESTCAL is similar to 
that adopted in CRKGRO, the detailed level crack-growth analysis program. 

A detailed description of the methodology is documented in the User's Manual 
of CRKGRO (5). The following paragraphs describe the highlights of the 
methodology. 


2.3. CRACK-GROWTH RATE EQUATION 
For cyclic loadings containing positive stress ratios (i.e., R 20), the 
modified Walker equation is used to describe the crack-prowth rate per cycle 


(da/dN}. The modified Walker equation can be expressed as: 


For AK > AK? R 20 


n 
a : >, 1-m = + ae 
da/dN = C [esa R) ,R< eee R 
R + R = + 
cut’? — cut 


For Ak $< AKy R 20 


da/dN = 0 


where C and n are the growth rate constants, m is the stress-ratio collapsing 
factor, and Rie is the cutoff values of positive stress ratio. 


For cyclic loadings containing negative stress ratios (i.e., R < 0), the 
Chang acceleration equat ion®) is used. In mathematical form, the Chang 
equation 1s expressed as: 


aD — 
da/dN = C|](1+R) k ,R2R R=R 


RR i ECR ae 


where q is the acceleration index, and Ro is the cutoff value for the nega- 


: : it 
tive stress ratio. 


2.4.0 LOAD INTERACTION MODEL 


- 


The Willenborg/Chang load interaction mode} | is adopted in ESTCAL to 
account for the spectnm loading effects. This model uses the generalized 
Willenborg mode | 8) to account for the tensile overload retardation effects, 
and the Chang acceleration scheme to account for the negative stress effects. 
The generalized Willenborg model can be written in the following form: 


1/2 
‘ Ol, Aa 
(kK Jie =k = Nov 5 aera - k 
max ett » max % » 
max ol. max 
Lhe 
ae oL Aa 
eee ae - ok ‘ee -K 
min eft we, max 7 w 
main ol max 
t= - (kK ) (Re 1) 
Max ee SO 
mech max : 
where ky. is the stress-intensity-factor corresponding to the maximum 


remotely applied stress, nol. is the stress-intensity-factor corresponding 
to the maximum stress of the overload, Aa is the incremental growth follow- 
ing the overload, Zoj, is the overload interaction zone size, and Reg is the 
overload shutoff ratio. 


For spectrum loadings, the effective stress-intensity-factor range and 
effective stress ratio are expressed in terms of the maximu., and minimum 
effective stress intensity factors as follows: 


AK Set ~ UN ote 7 (hninleff = Ako 
Ree ; CNien arte nak eee 


In the load-interaction-accounted-for option, the program uses the follow- 
ing equation to account for tensile overload retardation effect: 


- LAR = s: 
Eo} hate Rey Bor : 
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a ta a a, 


da/dN = 0 


where C, n, m, and Reat are the same crack-growth rate parameters described 
under "Fatigue Crack-Growth-Rate Equation.'" The threshold values of the 
stress-intensity-factor range are also identical to those used in the 
constant-amplitude cases. 


If the effective stress ratio is negative (i.e., Reff < 0), the Chang 
negative stress ratio equation is used in this program, which accounts for 
the compressive load acceleration effect: 


Z| 
t 


n 
5 
a = Do q & PR > = 
oe c|(1 * Ret) (Sa 7 RF pete eff Rage 
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where q is the acceleration index determined from test data generated for a 
specific negative stress ratio R < 0 and its R = 0 counterpart. 


Reduction of the overload retardation effect caused by a compressive 
spike load immediately following the tensile overload is accounted for 
through an effective overload interaction zone concept proposed by Chang (1), 
The effective overload interaction zone is defined in terms of the negative 


effective stress ratio (Roce < 0) as: . 


Pad 
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where Coz, is the plastic zone size introduced by the tensile overload. 


The plane strain plastic zone size is used if the stress intensity factor 
at the maximum depth for a part-through crack is to be calculated. The plane 
stress plastic zone size is used at the length direction for TC's and PTC's. 
The plane stress and plane strain plastic zone sizes are: 
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where Fy, is the material tensile vield strength. 


2.5. DAMAGE ACCUMULATION SCHEME: 


The Vroman linear approximation method hus been incorporated into tis 
computer program us the damage accumulation scheme. The following para- 
graphs brietly describe the method. 


For a given load spectrum as schematically shown in Table 1, the Vroman 
damage accumulation scheine proceeds by considering a load step (i) and using 
Taxi and Oming to calculate (da/dN);. The value of (0. 01a) /(da/dn) ; is then 
compared to Nj, where "a" is the instantaneous crack size. If (0.0la)/(da/dNt; 
is greater than Nj, then the crack growth for that particular load step is 
Oa = Nj x (da/dN);, the crack has then grown from "a" to (a + 4a), and the 
program proceeds to the next load step. 


If (0.01a)/(da/dN)j is less than or equal to Nj, the crack size will be 
fa + 0.0la), and this load step is reexamined. This process continues with 
(0.0la)/(da/dN)j being compared to the remaining cycles in the step. hen all 
load steps in the block or flight have been examined, the program then pro- 
ceeds to the first step of the next block (or flight) and continues. 


TABLE 1. A TYPICAL STRESS SPECTRUM TABLE (SCHEMATIC) 


Step Max Stress Min Stress 
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2.0. CRACK GROWTH UPDATING SCHEME 


The FLTGRO program uses a randem «ccess disk file to store crack-vrowth 
data for each control points and each type of mission. Hach control point 
requires the following information: (1) material constants, type of crack, 
crack geometry, and other parameters necessary to grow a crack, (2) mission 
type and associated prowth rate constants, and (3) a base mission mix. for 
the crack-growth-per-flight (da/df) approach, ys growth-rate parameter stored 
on the tracking operating file are C,a, wo-)! “, the stress ratio, and 
the maximum spectrum stress. The multiseyment-per-flight method stores vrowth 
rate parameter for each flight segment: i.e., c, n, the RMS of the maximum 
stress ina flight segment, the RMS ofopjq in a flight segment, and the nun- 
ber of Cycles per segment. The base mission mix is a unitblock of m missions 
and the number of occurrences per mission. 


In any time frame of the Individual Aircraft Tracking (AT) program, 
crack growth at a control point can be tracked by executing the FLTGRO pro- 
gram for the missions that have been completed and accessing the data base to 
obtain the growth rate parameters for each mission. Capability is provided 
for adding new missions to the data base or to the baseline mission mix. The 
crack-growth equation coefficients are computed and optionally replaced or 
appended to the data hase. If the missions flown were deviating from the 
baseline mission mix, the base mix can be updated to reflect the chanye. 


The FATGRO program produces a plot of crack length versus life in flights 
for the actual nu-ber of missions completed. The baseline mission mix can 
be optionally plotted on the same graph with the actual crack-growth history 
or can be displayed individually. The baseline mission is grown to failure to 
illustrate the life of the control point specimen. Another option is avail- 
able which appends the base mission crack-growth history the actual yrowth 
history. A combination of all these options in figure 2 shows the life of a 
control point it the base mission was flown, in comparison to the actual lite 
and the predicted life if the base mission is appended. 


CRACK LENGTH CINCHES) 


JATGRO TRACKING PROGRAM CRACK GROWTH ANALYSIS 
M-9) WITH DIFFERENT MISSION PARAMETERS 


FLTGRO Updating Crack Growth Analysis Results 


Section [1] 


PROGRAM OUTLINES 


The program identification is FLTGRO, FLTGRO consists of 14 routines. 
The functions of the main routine and 13 subroutines are: 


l. 


FLTGRO 


BFBD 


CENTER 


CLAMDA 


CYCCNT 


ESTCAL 


GROFLT 


LIST 


MANAGE 


NEWLEN 


PARAM 


PLOT 


TRP2 


WRALEC 


Supervisory routine, handles flow, checks keywords, gives 
diagnostics. 


Stores keywords for exact code titles and bach-face correc 
tion tables. 


Centers title cards for output. 


Computes and prints out C and + values for the da/dl rate 
equation in Method I, and C and 4 values for the da/ds 
rate equation in Method I]. 


Range-pair-counts the original spectrum loading. 


Reads input, prints output, prepares k-equations, yrons 
the crack using crack-growth rate per cycle (da/dN) 
equation for a series of initial crack size, aj, to the 
final size, ag. It also prepares tables of Cj, ‘j for 
mission-mix option, : ; 


Prepares K-equations, grows the crack using the crack- 
growth rate per flight (da/dF) equation from original to 
critical crack size. Prints out table of crack size versus 
flights. 


Prints out FLT records stored on the database. 


A data base managing routine to read and write FLT 
records. 


Brings user's file mmber on input to be recognized. 
Sets up parameters of analysis results for plotting. 
Plots grids, labels, and curves tor PARI. 


Parabolic interpolation routine for bach-face correst ion 
factor. 


Prints out crack-growth data for the series of 


Pilea, Sh -C.. srowths . 
initials “finals STO" 
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Section IV 


FLTGRO INPUT DATA DFCK 


The input data deck for FLTGRO is described in this section. A brief 
description of each type of input data card is presented in the following. 
The overall deck setup is shown in Figure 2, and a detailed description of 
each input card is given on the following pages. 


Card 1 FLT database disposition 


Card 2a ''Compute C and Lambda'' keyword 
2b Control point identifier 
2c Problem identification (i.e., title card) 
2d Material identification and properties 
2e Root-mean-square power exponent and shutoff ratio 
2f Stress ratio cutoff and crack-growth rate equation coefficients 
2g Initial crack sizes and geometry 
2h Crack coding control 
2i Number of initial crack sizes for computing C and A 
2j Initial crack sizes 
2k Method-type controls 
21 Mission identification 
2m Title identification for a given mission 
2n Limit stress and control parameters 
20 Mission loading and printing control parameters 
2p Stress spectrum 


Card 3a "Grow Crack'' keyword 

3b Control point identifier 
3c Growing crack control parameters 
3d Mission sequence definition 


Card da "Plot"! keyword 
4b Plotting control parameters 


Card SS "List" keyword 


Card 6 "End" keyword 
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CARD 6 
"END'' KEYWORD 


CARD 5 
ahrey” 


KEYWORD 


CARDS 4A-4B 
"PLOT'' KEYWORD 


CARDS 3A-3D 
"'GROW CRACK'' KEYWORD 
. eT RY 


-2P 
"COMPUTE C AND A'' KEYWORD 


CARD | 
FLT DATABASE DISPOSITION 


Figure 3. FLTGRO Input Deck Setup 


INPUT DATA CARD 1 


Description: Database disposition 


Format and example: 


Field Contents 
Option "FLTOF" keyword, beginning in column l. 
Disp Disposition of FLT database beginning in column 11. 


NONE - No FLT operating file will be used. 
NEW - A new file is to be created. 
OLD - An existing file is attached. 
Remarks 1. This must be the first data card in the data input deck. 


2. The Disp field is left-adjusted. 
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1 
{ 
| 
INPUT DATA CARD 2a 
Nesecription: Compute keyword card 
Format and exumple: 


1 20 St 
OPTION 
COMPUTE C AND LAMBDA 


Field Contents 
Option "Compute C and Lambda'' keyword, beginning in column 1. 
Remarks 1, This card initiates the input for the data needed to compute 
C anda. 


2. If C and A were previously computed and stored, an "old" 
file is attached; then, input cards 2a-2r will ke omitted. 


INPUT DATA CARD 2b 


Description: Control point identifier 


Format and example: 


] 8 18 80 


Field Contents 


CPIDNT Any alphanumeric information, maximum of eight characters, 
the user desires to input for identification of the data 
calculated by the "Compute C and lambda" option. 


LREP For old control points saved on the FLT file. 


BLANK - All requir 4 input data; cards 2c-2h will be read 
from the FLT file. 


REPLACE - A new set of input data, read from cards 2c-2h, 
will be read and stored on the FLT file replacing 
the old set of data. Left-adjusted, starting in 
column 9. 


For new control points, the input data, cards 2c-2h, will 


have to be provided in the input stream and IREP shall be 
left blank. 
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INPUT DATA CARD 2c 


Description: Problem identitication title card 
Format and example: 


Column 1 ite 


PROBID (i) 


SINGLE CYCLE PER FLIGHT, METHOD I, 30 kSI 


Field Contents 
PROBID (1) Any alphanumeric information which the user desires to input 


for problem identification. 


Remarks Only the first 72 columns are used for input. 


2] 


\ 
! 
i 
I 
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INPUT DATA CARD 2d 


Description: Material identification 
Format and example: 


Column 1 28 36 48 60 72 80 


Field Contents 
GEMATD Any alphanumeric information which the user desires to input for 
material identification (columns 1-36 only). 


BISLP Option to use one- or two-slope rate curve Blank - analysis 
with one-slope rate curve. '"'BISLOPE" - analysis with two-s!ope 
rate curve. 


ST GMAY Material vield strength (ksi) (F12.0) 


AKIC Plane strain fracture toughness, (ksi Vin.), used for part-through 
crack instability criterion on crack depth a. (F12.0) 


CKC Plane stress fracture toughness, (ksi Vin.), used for through- 
crack instability criterion on half-crack length c. (F12.0) 
Remarks : BISLOPE is a kevword, starting in column 29, 
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INPUT DATA CARD 2e 


Description: Crack growth rate equation coefficients, part 1 


Format and example: 


Column 


Field Contents 


POWER Power iran b to calculate the average of the b-th power 
product, (402)1/b. where ao is the stress range. (FLU.0) 


RETDA Retardation shutoff ratio. (F10.0) 


INPUT DATA CARD 2f 


Description: Stress ratio cutoff values and crack-growth rate equation 
parameters 


Format and example: 


Column 1 6 12 18 24 
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Field Contents 
RCUT The cutoff value of the positive stress ratio ''-R'' below which it 


it is assumed the material does not show stress ratio layering 
in (da/dN) vs AK plots. (F6.0) 


RCUTN The cutoff value of the negative stress ratio "'-R" below which 
it is assumed the material does not show further accelerated 
growth rate. (F6.0) 


EXPM Walker stress ratio effect exponent. (F6.0) 

EXPQ Chang growth-rate equation (for R < 0) exponent. (Fo.0) 
CWALK Walker growth-rate equation (for R 2 0) coefficient. (H12.0) 
EXPN Walker growth-rate equation (for R = 0) exponent. (F6.0) 
THA Threshold constant (F6.0) 


1. - Variable threshold 
0. - Fixed threshold 


DELKTH The threshold value of AK (ksi Vin.) obtained from R = 0 test. { 
(F12.0) 


eee 


INPUT DATA CARD 2f (Cone Luded ) 


Remarks : If "BISLOPE" option was used on card 2d, then an extra input 
card follows for the lower region of the rate curve: 


Extra input card for "BISLOPE" analysis 


a 12 _d8 3 
=XPQL Lavaca. 
iI .0 255-15) 12.9 | 
EXPMI. 
EXPOL As defined before, but for the lower part i.e.: Region I of 
CWALKI. the da/dn vs AK rate curve. 
EXPNL 
TRANSI. The AK value of the transition from upper to lower curve, 1.¢.: 
from Region I? to Region [. (F12.0) 
ATLEV The level of da/dn for the transition from Region II to Region 1, 


(E12.0) 


INPUT DATA CARD 2g 


Description: Initial crack sizes and geometry 
Format and example: 


Column 1 12 24 36 48 60 80 


Field Contents 
AINIT Initial value of crack depth a (inches) for part-through crack. 
(F12.0) 
ASPR Aspect ratio, (a/2c), for part-through crack. (F12.0) 
RADIUS Radius of hole (inch), if no hole leave blank. (F12.0) 
CINIT Initial value of crack length, c (inches); half of surface 


length dimension for centered cracks, full length for edge 
cracks. (F12.0) 


THK Equivalent thickness for transition; usually, it is the thick- 
ness of the plate. (F12.0) 
WIDTH Width of the plate, (inches). For centered cracks, the program 


’ 


halves the value of the width of the plate. (F12.0) 


INPUT DATA CARD 21 


Description: Number of initial crack sizes for computing C and aA 


Format and exanple: 


Column 
Field Contents 
NCASE Number of initial crack sizes to be read from card 2j. Nease 
is required for Methods I and I] and is right-adjusted. (15) 


INPUT DATA CARD 2h 


Description: Crack code controls 
Format and example: 


Column 1 6 12 80 


Contents 


Crack code no. (See crack library, p. ) 
Starting in colum 1 


Retardation option: (I6} 
0 = without load interaction 
1 = with load interaction 


Remarks Crack code no. left adjusted 
NRETRD, right adjusted 


INPUT DATA CARD 2) 


Description: Initial crack sizes 
Format and example: 


Column 1 5 10 1S 20) 
Fleer ope a ce es 


Field Contents 


The initial crack lengths, NC\SE lengths are read in. For part- 
through-cracks Ci's are the initial crack depths. For through- 
cracks Ci's are half of the initial surface lengths for centered 


> 


cracks; full lengths for edge cracks. (Refer to input curds 2¢ 
and 2h.) (FS.0) 


J 
is 


INPUT DATA CARD 2m 


Description: Method-type controls 
Format and example: 


Column 1 5 15 80 


Field Contents 


NPROB The number of different cases for which parameter c and n needed 
be estimated (15). 


METHOD 1 - Determine C and A for each mission (cycle/flight). 
2 - Determine C and A for each mission segment. (15) 


NSEG Number of segments in a flight to compute C and d for 
Method II only. (15) 


Remarks All input fields for this card are right-adjusted. 
1. Default is Method I. 
2. If NPROB >1, then NPROB sets of cards 2n-2y will follow 
when spectra are stored on files, or NPROB scts of cards 


2n-2r follow if spectra are on cards. 


3. A total of 50 sets of C and A can be determined. 
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INPUT DATA CARD 2n 


Description: Mission segment identification for loading condition 
Format and example: 


Column 1 8 


12 20 
MISID(i) |ISEG(i) | MISID (i+1) 
aasoxst [| 


24 
jisr6 (+1) | 
Heeaae © 


80 
Ion fa Iwisw ose [ise oe [id 
ek ee ee erates vee ae 


Field Contents 


For Method I: 


MISID (i) Mission identifier (maximum of eight characters); rest of 
¢ard is blank). (A8) 


For Method II: 
MISID (i) Mission segment identification (maximum of eight characters). 


ISEG (i) Terminating cycle number in flight for mission segment, right- 
adjusted. (14) 


ISEG (NSEG) is equal to the total number of cycles in flight. 


NSEG is given on card 2m, 


INPUT DATA CARD 2o 


Description: Title identification 


Format and example: 


Column 1 68 80 | 
Field Contents 
NAME Title identification (up to 68 characters) for mission defined 


by cards 2m and 2n. 
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INPUT DATA CARD 2p 


hescription: Lunit stress and control parameters 
Format and example: 


SU 


Column l 12 1? 29 
SIGLIM MOSGMA } SMULTI 
Ecco Ee | [ 


Field Contents 
SIGLIM Limit stress (ksi). (F12.0) 
MOSGMA Stress input control parameter (I5) 


20 - stress is in terms of percent of limit stress. 
Q - stress is in terms of ksi. 


SMULTI Stress magnification or reduction factor to bring spectrum to 
desired stress level.  (F12.0) 


Remarks Mosema is right-adjusted. 
MOSOMA 2 0 Stress = SSE S Smith 
iN MOSGMA < 0 Stress = Stress x Smulti 
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INPUT DATA CARD 2q 


Description: Mzssion loading and printing control parameters 


Format and example: 


Column 
Field Contents 
NFLIGHT Number of unit blocks applied for crack growth computation to — 
estimate C, A and product. 
L Number of cycles in a unit-block, or given subunit-block. 
NUNIT Number of flights in a unit-block or given subunit-block. 
MOPEAK Input peaks-valleys adjustment parameter: 
20 - MIN-MAX input 
< 0 - MAX-MIN input 
MOCRK Print control parameter for crack growth data: 
20 - Crack-growth data are printed. 
< 0 - Crack-growth data are not printed. 
KCON Print control parameter for stress history: 
> 0 - Input stress history is not printed. 
s 0 - Input stress history is printed. 
NFILE Alternate file where stress spectrum is stored. This file 
input unit number must be between 10 and 98, inclusive. If 
NFILE is negative, this unit is rewound before reading of 
file. If 0 or 5, card input is expected. 
NRP Control for range-pair counting: 
0 = No range-pair counting. 
1 = Range-pair counting. 
Remarks All input fields for this card are right-adjusted. 


INPUT DATA CARD 2r 


Description: Stress spectrum ' 


Format and example: 


Column 1 5 15 25 35 
DATA sini fovciei [2 
ee Oe, 
or 
1 5 15 25 35 | 
DATA | SMINi SMAXi CYCLF1 | | 
Pfs isn io 3 
Field Contents 
DATA Blank on all SMAXi, SMINi, CYCLEi cards. 
'"End'' on extra card marking the end of a spectrum, starting i 
in column l. 
SMIN Minimum stress value (ksi). 
F SMAX Maximum stress value (ksi). 
Remarks 1. The spectrum is read on the unit defined by NFILF on 


card 2q. 


2. The min-max or max-min input is defined by MOPEAK on 
card 2q. 


remem ee = 
~) 


uw 


For mission mix option, spectrum is limited to 3000 or less 
steps (i) in each uniblock. 


INPUT DATA CARD 3a 


Description: "Grow Crack"' keyword card 


Format and example: 


80) 


Column 1 10 


Field Contents 
OPTION "GROW CRACK'' keyword, beginning in column 1. 


{ 
i 
\" 


= . 
2a ene cea creat ti 


i, ie 


INPUT DATA CARD 3b 


Description: Control point identifier 


Format and example: 


Column 1 8 80 
SU sa ete nee et 
M-85 
Field Contents | 
CPIDNT The alphanumerical information (maximum of eight characters) 


used for identification of the data that were stored is now 
attached and is to be read from the FLT data base. This card 
is required only when the "compute C and lambda'' option is not 
used. 


INPUT DATA CARD 3c 


Nescription: Growing crack control parameters 
Format and example: 


Column 1 8 16 


Ww 
le) 


Field. Contents 
CF Final crack length (inches). 
NBLKS Number of blocks (limiting number of flights for the C vs flight | 
growth curve). 
MMIX Number of mission strings. 
NBASE =] - If final or critical crack size is not reached upon 


completion of the discrete missions mix, append the 
base mission mix. 
= 0 - Only discrete mission mix is used. 


NREPET = ] - Grow the crack using the base mission mix. 
= 0 - No operation. 


NBASEI = ] - [Input mission mix is to he stored as the base mission mix, 
= 0 - No operation. 


NRETRD = 1 - Load interaction is considered (for Method II only). 
= 0 - Load interaction is not considered. 


FNTBLK Total number of flights in one mission (for Method II only). 


Remarks All input fields for this card are right-adjusted, except for 
CF. 


Description: 


INPUT DATA CARD 3d 


Mission sequence definition 


Format and example: 


Column 


Field 
NFAC 


MIS 


Remarks 


1 8 12 


NFAC(i) [MIS(i) a asian _[oene 


ss NFAC(it+3) | MIS(i NFAC(i+4) ca 


1 
56 60 68 72 
S [ars (ira) NFAC (i+5) MIS(it+S) 


Ww 
ty 


Pee 


80 


Contents 


The number of times the loading condition is to be repeated. 
The sequence number for the loading condition as saved in the 
FLT database. The sequence number can be obtained from the 
printout of the "LIST" option. 

1. All input fields for this card are right-adjusted. 


2. Card 3d is needed for Methods I and II. 


3. Mission mix applies to Method I only. 


INPUT DATA CARD 4a 


Description: Plot keyword card 
Format and example: 


Column 1 80 


PLOT 
Field Contents 
OPTION "PLOT keyword, beginning in column 1. 
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Description: 


INPUT DATA CARD 4b 


Plotting control parameters 


Format and example: 


Column 


Field 


OPT 


SCALE(2,1) 


Remarks 


1 5 10 15 20 25 
SCALF(1,1) |S ; 
| 
eee 
25 30 80 
Contents 


Specifying parameters to be plotted. 

1 - Plot 

0 - No plot 

OPT = Crack size versus life in flights. 


Array specifying linear, semilog, or log-log grid scaling. 
SCALE (1,) = 0 - x-axis is linear 

1 - x-axis is log 
SCALE (2,) = 0 - v-axis is linear 

1 - y-axis is log 


All values are integers, right-adjusted. 


INPUT DATA CARD 5 


Description: LIST keyword card 
Format and example: 


Column 1 80 


Nn rr ee 


Field Contents 
OPTION "LIST" keyword, beginning in column 1. 
Remarks This option will list the loading conditions that have been 


saved on the FLT database. 


INPUT DATA CARD 6 


Description: END keyword card 
Format and example: 


Column 


Field Contents 


OPTION "END" keyword, beginning in column 1. 


Remarks This card terminates the input data deck and signals the 
program to terminate execution. 


- 


Section V 
EXAMPLE CASES 


This section presents two example cases which are designed to illustrate 
the capability of the FLTGRO program. The first example is to predict the 
crack-growth behavior of a 2219-T851 aluminum center-craked-tension (CCT) 
specimen subjected to a typical figher spectrum loading. Figure 4 shows the 
specimen dimensions and the crack configuration by the single-cycle-per-flight 
method (Method I). 


The crack-growth rate constants and parameters used in the example are 


as follows: 
_ -10 ae 
C = 5.066 x 10 Rog 3 
n= 3,83 q= 1.0 
: = +0.75 
m = 0.6 cae ; 
ann te - ~~ = -0.99 
Okey 1.5 ksi. in. ce 
ae 43 = 48 ksi 
Ke 65 ksi in. Cry 
bs 2 


The spectrum used in this example was a typical fighter aircraft air-to- 
ground (A-G) mission, generated in phase III of this research and development 
effort (4). The detailed spectrum table is shown in the appendix. The spectrum 
is in the random cycle-by-cycle format which contains 150 flights, each flight 
consisting of 19 cycles. Each peak and valley is in the percentage of design 
limit stress (% of DLS). The random spectrum was range-pair counted before it 
was used in the analysis. The range-pair counted spectrum was printed out as 
shown in the example. All the peaks and valleys were converted to stress in 
KSI unit, based on o,._ = 24 ksi. 

lim { 

Eight (8) initial crack sizes were selected in order to characterize the 

A-G mission. These 8 initial crack sizes were: 


C, = 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.9, and 1.1 (all in inches). 


The input echoes are shown in the following computer print-out. A table 
of dc/dF vs K for these 8 points was also prepared by the FLTGRO program 
shown in the output printout. The characterized A-G mission crack-growth- 
rate-per-flight parameters obtained by FLTGRO are: 


C = 2,0648 x 10° 


A 


3.6119 
= 10.692 ksi 


These characterized crack-growth-rate-per-flight parameters were then 
used to perform crack-growth prediction on a case. The initial crack size 
was Cj = 0.145 inch. All crack growth parameters used in the life prediction 
are listed in the printout. The crack-growth summary table is shown in the 
last page of the computer printed outputs of this example. The table shows 
the FLTGRO predicted life from C; = 0.145 in. to failure was Np = 5462 flights. 


Section A-A 


ae 


gas 2G; 


Material 
2219-T851 aluminum plate 


(—— ne 6 ee 


Figure 4. Test Specimen Configuration 
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The second example is to predict the crack-growth behavior of the 

rT) identical crack configuration under the identical spectrum loading as that 
was described in the first example, by the multi-segment per-flight method 
(Method II). The crack growth rate constants and parameters used in the 
second example are identical to those used in the first example. Each flight 
of the A-G spectrum was arbitrarily divided into three segments, the first 
segment, (A-G), consisting 5 cycles, the second segment, (A-G)7 consisting 
7 cycles, and (A-G)z also consisting of 7 cycles. 


Again, eight (8) initial crack sizes were selected for the character- 
ization of each flight segment. The initial crack sizes selected for this 
example were identical to those in the first example. The calculated C.d 
and @ for each flight segment are: 


(A-G)]: C = 4.072 x 1079, d= 3.7664, a= 10.692 Ksi 
(A-G)o: C= 5.2411 x 1079, d= 3.7592, «= 10.692 Ksi 
(A-G)z C= 6.5638 x 1079, h= 3.7634, o= 10.692 Ksi 


These parameters were then used in the prediction of the crack growth 
of a CCT specimen with an initial size C = 0.145 in. The crack growth 
results are shown in the summary table of the computer outputs. It shows 
the predicted life by Method II is Np = 6246 flights. A CRT plot of the ‘ 
crack growth curve is shown in the printout. 
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Appendix 
A TYPICAL FIGHTER AIR-TO- GROUND 
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